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ABSTRACT 


Detailed  measurements  were  performed  in  the  region  of  interaction 
of  a  laminar  boundary  layer  with  a  compression  corner  at  Mach  numbers 
near  2.5.  Different  models  were  tested  permitting  variation  in  the 
angle  of  compression  and  the  conditions  downstream  of  reattachment. 

The  heat  flux  and  the  resistance  of  equilibrium  of  a  hot  wire 
anemometer  at  any  location  in  the  flow  field  were  measured. 

These  measurements  were  supplemented  by  either  the  wall  or  the  pitot 
pressure  to  compute  all  thermodynamic  and  dynamic  variables. 

The  results  permitted  verification  of  the  hypothesis  of  zero 
normal  pressure  gradients  when  the  compression  angle  is  small  and  the 
boundary  layer  is  laminar.  The  enthalpy  in  the  separated  bubble  was 
constant  and  equal  to  the  wall  enthalpy  of  a  flat  plate  with  attached 
laminar  boundary  layer  at  the  same  Mach  number.  Evidence  of  a  reversed 
flow  with  velocities  approximately  5%  of  the  free  stream  velocity  was 
observed.  The  recompression  along  the  separating  streamline  was  found 
to  be  very  nearly  isentropic.  Critical  points  were  located  in  the 
reattachment  region  and  a  physical  explanation  proposed  with  special 
emphasis  on  the  location  of  transition. 
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1.0  INTRODUCTION 

For  the  design  of  present  day  high  speed  aircraft  it  is  becoming 
increasingly  important  to  understand  the  behavior  of  supersonic  boindary 
layers  when  exposed  to  sudden  changes  in  free  stream  conditions.  These 
changes  could  be  caused,  for  example,  by  a  shock  impinging  on  the  boundary 
layer  or  by  a  rapid  change  in  the  body  shape.  Experience  has  shown  that 
even  though  the  region  where  the  viscosity  acts  is  very  thin,  under  cer¬ 
tain  flow  condi cions  this  layer  can  move  away  from  the  body,  creating 
large  shifts  of  the  aerodynamic  loads.  It  is,  therefore,  important  to 
be  able  to  predict  these  flow  conditions  and  their  effects  on  the  boundary 
layer. 

This  type  of  phenomenon  occurs  on  wings  with  trailing  edges  of  finite 
thickness,  causing  an  additional  drag  (base  drag);  it  also  occurs  on  re¬ 
entry  vehicles  with  flares,  introducing  a  favorable  effect  on  the  aero¬ 
dynamic  heating.  The  type  of  flow  investigated  in  this  paper  is  the 
separation  of  a  boundary  layer  caused  by  a  sudden  compression  such  as 
that  on  wings  with  flaps  or  ailerons,  or  the  conical  flare  which  follows 
the  cylindrical  section  of  a  Gemini  capsule.  These  flows  all  belong  to 
a  class  of  problems  called  separated  flows.  The  problem  which  interests 
us  is  one  of  the  most  difficult,  for,  as  will  be  seen  later,  neither  the 
point  of  separation  nor  that  of  reattachment  is  fixed. 

The  physical  features  of  the  flow  are  as  follows.  A  supersonic  stream 
flowing  past  a  flat  plate  is  suddenly  turned  through  a  corner  at  a  certain 
point  0  (see  Fig.  1).  If  the  effects  of  viscosity  are  neglected  the 
problem  is  a  trivial  one.  A  shock  originating  at  0  will  turn  the  stream 
through  the  deflection  angle  a.  If,  however,  viscosity  is  taken  into 
account  we  know  that  there  exists  a  thin  layer  near  the  wall  which  is 
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slowed  down  by  friction.  The  bottom  part  of  this  layer  will  be  subsonic, 
hence  a  shock  cannot  exist  there.  The  pressure  jump  will  be  spread  over 
a  finite  length.  When  this  layer  approaches  the  corner  It  will  sense  the 
pressure  Increase  before  reaching  0.  If  the  momentum  of  the  particles 
In  the  boundary  layer  Is  not  sufficiently  large  to  overcome  this  unfavor¬ 
able  pressure  Increase  the  particles  will  move  away  from  the  wall.  There 
then  exists  on  the  wall  a  stagnation  point  (S),  which  Is  called  the 
separation  point.  A  streamline  (SA)  originates  from  this  point.  As 
the  flow  separates  the  particles  will  pick  up  sufficient  momentum  to  bring 
them  into  contact  with  the  body  at  a  point  R,  which  Is  called  the  re¬ 
attachment  point.  This  point  Is  also  a  stagnation  point,  and  a  streamline 
A'R  terminates  there.  Using  simple  arguments  based  on  conservation  of 
mass  we  see  that,  provided  no  flow  crosses  the  solid  boundary  between  S 
and  R,  the  streamline  originating  in  S  Is  the  same  one  that  terminates 
at  R.  Hence  there  exists  a  closed  region  SAA'RO  which  always  consists 
c*  the  same  particles.  Looking  closely  at  this  region  we  see  that  the 
forces  acting  on  the  particles  within  It  are:  shear  stresses  on  the 
Interface  SR,  and  wall  friction  on  SRO.  The  particles  near  SR  will 
be  moving  downstream,  so  that  other  particles  will  have  to  replace  them. 
This  region  therefore  has  a  clockwise  circulating  motion.  Actually,  It 
Is  not  correct  to  call  the  line  SR  an  Interface,  because  a  priori  no 
discontinuities  of  any  sort  occur  on  this  line;  It  is  merely  the  line 
that  separates  the  particles  which  have  sufficient  momentum  to  overcome 
the  pressure  rise  at  the  reattachment  region  from  those  which  have  In¬ 
sufficient  momentum  to  overcome  the  pressure  rise  and  are  turned  back 
Into  the  region  SRO;  for  this  reason  SR  Is  called  the  separating 
streamline.  The  above  slmpl*  arguments  give  a  good  picture  of  the 
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expected  flow;  howe  ,  the  difficulty  lies  In  determining  the  locations 
of  S  and  R,  In  being  able  to  predict  this  location,  In  finding  SR, 
and  more  generally,  In  being  able  to  predict  what  will  happen  under  a 
given  set  of  conditions. 

Theoretically,  the  difficulties  appear  quite  discouraging  at  first 
glance.  The  principal  difficulty  Is  that  the  boundary  layer  equations 
are  not  expected  to  hold  at  the  point  of  separation,  since  a  normal 
pressure  gradient  may  exist  there.  Even  by  assuming  that  there  Is  In 
fact  no  appreciable  normal  pressure  gradient  It  Is  still  Impossible  to 
solve  these  equations  because  the  streamwise  pressure  distribution  Is 
one  of  the  unknowns.  However,  a  great  contribution  has  been  made  In  this 
regard  by  Chapman,  et  al . , 1  who  suggested  a  way  to  overcome  this  diffi¬ 
culty  by  assuming  that  the  pressure  distribution  is  conditioned  by  the 
interaction  between  the  laminar  boundary  layer  and  the  supersonic  free 
stream.  This  assumption  Is  justified  by  several  experiments  and  Is  now 
widely  used.  It  Is  called  the  "free  Interaction  hypothesis."  More  pre¬ 
cisely,  this  hypothesis  Implies  that  the  interaction  between  the  viscous 
layer  and  the  supersonic  free  stream  depends  solely  on  local  quantities 
and  Is  not  directly  dependent  upon  the  geometry  which  Induces  separation. 

An  extension  of  this  assumption  for  the  nonadiabatlc  wall  has  been 
2 

proved  by  Curie. 

Several  theoretical  solutions  were  found  by  using  this  simplifica¬ 
tion.  An  approximate  solution  based  on  the  application  of  Howarth’s 
transformation  to  the  compressible  boundary  layer  equations  has  been 
found  by  Curie.  Other  approximate  methods  based  on  consideration  of 
the  boundary  layer  as  two  layers,  an  Inner  layer  in  which  viscosity  and 
pressure  forces  are  dominant,  and  an  outer  layer  with  Inertia  forces, 
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were  also  found  by  Gadd,  and  Hakklnen,  et  al.  More  sophisticated 

S  6 

mathematical  methods  were  used  by  Holt  and  Lees  and  Reeves.  Holt  uses 
a  multi -moment  method  Instead  of  edge  boundary  conditions  as  in 
Pohlhausen's  method  and  expresses  the  shear  stress  as  a  polynomial  In 
u  as  opposed  to  using  a  velocity  polynomial.  The  method  of  Lees  and 
Reeves  uses  a  combination  of  locally  similar  and  momentum  techniques. 

Both  solutions  use  a  coupling  equation  between  the  dissipative  flow  and 
the  Isentropic  free  stream  to  account  for  the  free  interaction.  Holt's 
method  of  solution,  called  the  "method  of  integral  relations,"  has  the 
advantage  of  simplicity. 

Experimentally  this  problem  has  also  attracted  the  attention  of 
several  researchers.  One  of  the  earliest  and  still  very  important  con¬ 
tributions  Is  the  work  of  Chapman,  Kuehn  and  Larson1  in  which  several 
other  separated  flows  are  examined.  In  this  work  experimental  evidence 
of  the  free  interaction  hypothesis  is  presented,  and  the  influence  of 
Mach  and  Reynolds  numbers  on  separated  laminar  transitional  and  turbulent 
boundary  layers  is  shown  in  great  detail.  Several  other  investigations 

examining  different  aspects  of  the  problem  were  carried  out  by  Gadd,^ 
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Bogdonoff  and  Vas,  Needham.  Recently  Lewis  extended  Chapman's  free 
Interaction  correlation  to  hypersonic  adiabatic  and  cooled  walls  using 
Curie's  approach.  Although  a  slightly  different  problem  is  treated  by 
Sirlcx,  Mirande  and  Delery,11  we  will  mention  their  work  which  is  con¬ 
cerned  with  a  detailed  Investigation  of  the  reattachment  region.  To 
Illustrate  the  complexity  of  this  type  of  flow  we  will  also  mention 
Glnoux's  paper  which  shows  the  existence  of  three  dimensional  effects 
(generally  Interpreted  as  Goertler's  vortices)  In  flows  which  are  essen¬ 
tially  two-dimensional. 
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With  the  impressive  amount  of  research  already  undertaken  in  this 
field  the  present  experimental  investigation  finds  its  justification  in 
the  fact  that  most,  if  not  all,  of  the  preceding  work  consisted  of 
examining  global  flow  quantities  with  little  attention  paid  to  detailed 
properties.  For  example,  to  our  knowledge  there  has  been  no  final 
experimental  evidence  of  the  existence  of  a  reversed  flow,  although 
this  is  physically  unquestionable.  We  also  feel  that  the  computational 
difficulties  encountered  by  most  theoreticians  in  extending  the  solution 
to  the  reattachment  region  may  not  be  purely  mathematical,  but  might  be 
due  to  the  fact  that  the  assumptions  used  are  not  valid  in  that  region. 

Since  our  purpose  was  to  provide  more  detailed  knowledge  of  separa¬ 
tion  phenomena,  we  sought  to  measure  the  thermodynamic  and  dynamic 
quantities  throughout  the  flow  field  using  hot  wire  and  pitot  probes. 
Since  the  recompression  region  is  one  of  the  least  understood  anas  we 
also  attempted  to  localize  the  critical  points  in  that  region  and  to 
interpret  their  roles  in  terms  of  local  quantities  of  the  boundary  layer 
in  that  region,  following  the  work  of  Siriex,  et  al J 

During  t  >sent  study  the  Mach  and  Reynolds  numbers  were  varied 
only  slightly  .  .  latter  always  having  a  value  low  enough  to  Insure 
laminar  flow  over  the  whole  interaction  region.  The  models  tested  were 
all  adiabatic  and  allowed  for  variation  of  the  compression  angle  a  and 
the  downstream  conditions.  The  flow  was  carefully  checked  for  two- 
dimensional  behavior. 


6 


2.0  EXPERIMENTAL  APPARATUS 

2.1  Wind  Tunnel 

The  experiments  reported  here  were  carried  out  in  the  Aeronaut¬ 
ical  Sciences  Division  6"  x  6"  supersonic  wind  tunnel  at  the  University 
of  California  in  Berkeley.  This  facility  (described  in  extensive  detail 
by  Bossel  )  is  a  closed  type  continuous  flow  tunnel.  Stagnation  temper¬ 
ature  and  pressure  are  adjustable  in  the  range  150  <  T^  <  50°F  and 
1.3  <  P^  <  35  psia,  which  allows  a  variation  of  the  Reynolds  number  per 
foot  in  the  range  1.50  x  10^  <  Re/ft  <  7.35  x  10^  at  a  Mach  number  of 
2.60.  The  free  stream  Mach  number  in  the  test  section  varies  by  3%  in 
the  streamwise  direction  over  a  distance  of  the  same  order  of  magnitude 
as  the  length  of  the  model. 

This  facility  is  equipped  with  a  schlieren  and  shadowgraph 
system  which  permits  ore  to  observe  and  photograph  the  flow. 

2.2  Models  Used 

Two  models  were  used  for  the  investigation,  both  of  which  were 

adiabatic. 

Model  A  (see  Fig.  2)  was  essentially  used  to  check  for  estab¬ 
lishment  ot  the  flow.  It  permits  continuous  variation  of  the  angle  a 
from  7°  to  23°.  It  also  allows  for  changes  in  the  distance  from  the 
leading  edge  to  the  compression  corner  from  3  to  4-1/2  inches,  and  could 
be  used  as  a  flat  plate  model  by  removing  the  ramp  mounting.  Pressure 
taps  at  1/8  inch  intervals  were  provided  on  this  model. 

Model  B  (see  Fig.  3)  had  the  same  provision  for  pressure  measure¬ 
ment;  however,  the  ramp  position  was  fixed  at  4  inches  from  the  leading 
edge.  The  compression  angle  could  be  set  at  the  values  9°,  11°,  13°,  and 
15°  by  changing  the  ramp  mounting.  This  model  also  had  an  additional 
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feature--the  ramp  was  constructed  of  a  flexible  copper  sheet  soft-soldered 
on  lateral  ribs  spaced  1/8  Inch  apart,  which  allowed  a  continuous  or  a 
sudden  deformation  at  any  angle  and  at  any  distance  from  the  compression 
corner.  This  deformation  could  be  produced  during  the  run.  This  feature 
was  used  to  enable  one  to  change  downstream  conditions  at  any  point  during 
and  after  reattachment. 

The  leading  edges  of  both  models  were  checked  carefully  and 
regularly  under  a  microscope;  this  proved  to  be  beneficial  In  establish¬ 
ing  a  laminar  boundary  layer. 

The  pressure  taps  on  both  models  were  connected  to  a  Decker 
differential  pressure  transducer  with  a  range  of  ±3  Inches  of  water, 
which  was  regularly  calibrated  against  a  Dibutyl -phthalate  micromanometer. 
The  voltage  output  of  the  transducer  system  was  recorded  on  a  Llbrascope 
X-Y  recorder. 

2 . 3  Probes 

Several  hot  wire  and  pitot  probes  were  tested  and  calibrated; 
however,  all  the  data  reproduced  in  this  report  were  obtained  with  the 
probes  shown  in  Fig.  4.  These  probes  were  chosen  because  they  Indicated 
a  minimum  interference  with  the  flow.  Both  probes  had  the  same  mounting 
--a  conical  section  terminated  by  a  cylinder  which  fits  the  arm  of  the 
traversing  mechanism. 

The  pitot  probe  consisted  of  a  piece  of  stainless  steel  tubing 
fitted  into  the  cone  at  an  angle  of  30°  and  bent  parallel  to  the  flow 
1/2  inch  from  its  tip.  The  portion  of  the  tube  at  an  angle  to  the  flow 
is  flattened  to  present  minimum  disturbance;  the  mouth  of  the  tube  is 
also  flattened  for  boundary  layer  measurements  (see  dimensions  in  Fig.  4). 
It  was  carefully  constructed  and  regularly  checked  under  a  microscope. 
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The  hot  wire  probe  had  a  thin  plate  fixed  to  the  cone  on  one 
er.d  which  held  one  needle  parallel  to  the  flow  on  each  side  at  the  other 
end.  This  plate  was  made  of  two  sheets  of  copper  glued  together  but 
Isolated  electrically  so  that  with  the  needles  soft  soldered  on  each 
face,  It  served  as  an  electric  conductor  between  the  hot  wire  and  the 
leads  to  the  electronic  equipment.  This  plate  could  be  made  very  rigid 
and  still  very  thin  (see  Fig.  4). 

Both  probes  were  mounted  on  the  test  section  traversing  mech¬ 
anism  which  permitted  positioning  on  the  three  axes  with  an  accuracy  of 
±0.001  Inch. 

The  pitot  probe  pressure  signal  was  connected  to  a  Decker  or  a 
Wlanko  ±5  psl  pressure  transducer  and  the  electric  signal  was  recorded  on 
the  X-Y  recorder  where  the  Y  axis  was  actuated  by  a  voltage  Input  pro¬ 
portional  to  the  distance  from  the  wall  of  the  probe  tip. 

The  hot  wire  was  connected  to  a  heat  flux  constant  temperature 
anemometer.  When  the  equipment  was  run  in  the  cold  resistance  mode  the 
resistance  of  the  unheated  wire  could  be  measured  with  an  accuracy  of 
better  than  ±0.01  fi.  Two  other  readings  from  the  wire  voltage  were  re¬ 
corded  for  very  small  overheats  and  since  the  bridge  was  at  balance  the 

_5 

wire  currents  could  be  computed  with  an  accuracy  of  better  than  ±3  x  10 


amps. 
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3.0  EXPERIMENTAL  TECHNIQUES  AND  DATA  REDUCTION 
3. 1  Flow  Establishment 

The  flow  to  be  investigated  Is  the  laminar  two-dimensional 
flow  over  a  compression  corner.  We  must  then  make  sure  the  two  condi¬ 
tions  are  satisfied,  namely:  two-dimensional  and  laminar  flow.  One 
additional  condition  is  that  the  boundary  layer  approaching  the  corner 
should  be  of  a  Blasius  type.  In  other  words,  since  our  model  has  a 
leading  edge  at  a  finite  distance  from  the  corner,  we  must  make  sure 
that  the  leading  edge  effect  has  died  out  before  the  interaction 
occurs. 

3.1.1  Two-dimensional  flow 

The  width  of  the  test  section  is  5-1/2  Inches,  which 

limits  the  aspect  ratio  of  the  model  particularly  because  the  length 

l  (see  Figs.  A  and  B)  cannot  be  made  shorter  than  a  certain  length 

JL.  for  the  reasons  explained  in  Section  3.1.2.  Due  to  these  llmita- 
mi  n  r 

tions  our  aspect  ratio  (s/ft)  was  of  the  order  of  1,  which  made  it 
important  to  investigate  the  two-dimensionality  of  the  flow.  An  ex¬ 
tensive  study  of  this  kind  has  been  performed  by  Lewis, ^  who  found 
that  by  mounting  side  plates  at  equal  distances  from  the  center  line 
and  by  varying  this  distance,  the  flow  tends  to  a  limit  when  the  aspect 
ratio  becomes  of  the  order  of  or  larger  than  one.  Using  different  wind 
tunnel  facilities  and  models,  Lewis  proved  that  this  limit  Is  actually 
the  desired  two-dimensional  flow.  This  result  is  used  in  the  present 
investigation.  However,  since  the  Mach  number  used  in  our  case  Is  of 
order  2.5,  an  aspect  ratio  study  was  carried  out.  When  AR  became 
equal  to  or  larger  than  one,  the  limit  was  attained  for  a  deflection 
angle  a  £  15°.  For  larger  a  a  considerably  larger  AR  was  needed. 
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As  proved  by  Lewis,  we  will  assume  that  the  limit  found  does  represent 
the  two-dimensional  flow. 

All  of  our  final  data  were  taken  on  model  B  where  a 
was  always  less  than  15°.  It  was  observed  (see  Fig.  5)  that  taking  out 
the  sideplates  altogether  does  not  substantially  affect  the  flow,  pro¬ 
vided  (s/A)  >  1.  The  difference  in  the  pressure  distribution  observed 
in  that  figure  near  the  end  of  the  ramp  is  due  to  the  appearance  of 
transition  when  the  sideplatcs  were  mounted  (see  subsequent  section). 

It  was  also  observed  that  transition  occurs  more  easily  when  the  model 
spans  the  wind  tunnel,  presumably  because  the  tunnel  wall  boundary  layer 
is  turbulent  and  this  turbulence  is  fed  to  the  main  flow  through  the 
separated  bubble  in  which  the  velocity  is  very  small. 

In  the  light  of  these  remarks,  it  was  decided  to  use 
model  B,  which  does  not  span  the  test  section  and  has  an  aspect  ratio 
of  1. 

A  final  check  was  performed  by  applying  a  thin  oil  film 
on  the  model  and  observing  the  oil  traces.  This  showed  that  there  exists 

a  region  roughly  3/4  to  1  inch  on  each  side  of  the  centerline  where  the 

oil  traces  were  essentially  parallel  to  the  main  stream.  Changing  the 
ramp  configuration  downstream  of  the  reattachment  location  did  change  the 
oil  traces,  but  the  region  of  two-dimensional  flow  was  still  present.  It 
will  be  assumed  in  what  follows  that  the  ramp  length,  or  more  precisely 
the  conditions  downstream  of  reattachment,  do  not  affect  the  two-dimen¬ 
sional  nature  of  the  flow.  This  is  an  important  point,  because  the  down¬ 
stream  condition  effect  on  the  flow  is  an  important  part  of  this  study. 

3.1.2  Boundary  layer  ahead  of  interaction 

As  stated  earlier  it  is  important  for  our  study  that  the 
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boundary  layer  ahead  of  the  Interaction  be  a  self-preserving  Blaslus 
flow.  Model  A,  which  has  provision  for  changing  the  distance  l  was 
used  first  without  a  ramp  and  then  with  the  ramp  at  different  locations 
l.  Pressure  distributions  In  the  streamwlse  direction  were  taken  (see 
Figs.  6),  and  it  was  found  that  l  *  4.00  inches  was  sufficient  to 
achieve  this  condition.  This  length  was  used  on  Model  B,  and  It  seems 
to  be  the  optimum  length  because  L  >  4.00  can  have  adverse  effects  on 
the  location  of  the  transition  (see  next  section). 

3.1.3  Laminar  flow 

The  last  condition  to  be  satisfied  is  that  the  flow  be 
laminar  before,  as  well  as  after,  the  interaction. 

It  is  a  known  property  of  laminar  boundary  layers  that, 
when  observed  by  a  shadowgraph,  they  exhibit  a  focusing  effect  that  pro¬ 
duces  a  bright  line  near  their  edge.  This  criterion  has  been  used  by 
Chapman  et  al J  in  their  investigation.  This  effect  is  generally  Inter¬ 
preted  by  the  fact  that  laminar  boundary  layers  have  a  density  profile 
of  small  curvature  near  the  edge  as  compared  with  a  turbulent  density 
profile  which  exhibits  a  fuller  curvature.  This  interpretation  was 
checked  by  hot  wire  measurements  on  a  flat  plate  boundary  layer.  Figure 
7  shows  two  temperature  profiles;  the  curvature  of  the  laminar  profile 
is  much  smaller  than  that  of  the  turbulent  one.  These  two  profiles 
exhibited  the  behavior  explained  earlier  when  observed  with  a  shadowgraph, 
and  a  confirmation  of  their  laminar  and  turbulent  natur.-  was  obtained  by 
a  hot  wire  fluctuation  qualitative  study.  These  temperature  profiles 
were  obtained  from  the  steady  state  hot  wire  measurements  explained  In  a 
subsequent  section. 
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Several  similar  tests  were  carried  out;  all  confirmed 
the  existence  of  a  bright  line  around  the  edge  of  a  laminar  boundary 
layer  when  observed  by  a  shadowgraph.  This  property  was  used  in  all 
subsequent  runs  where  we  could  decide  on  the  state  of  the  boundary  layer 
with  a  quick  look  at  the  shadowgraph  screen. 

In  order  that  transition  does  not  occur  in  the  region 
of  Interest  in  our  flow,  the  tunnel  was  used  at  its  minimum  stagnation 
pressure;  this  also  had  another  advantage,  because  under  these  condi¬ 
tions  the  stagnation  pressure  regulating  control  was  extremely  stable 
and  Identical  free  stream  conditions  could  be  reproduced  in  different 
runs.  When  transition  was  desired,  a  very  small  increase  of  stagnation 
pressure  would  move  the  transition  upstream  to  the  reattachment  region. 
The  stagnation  temperature  was  also  set  at  Its  minimum  value  for,  even 
though  this  had  an  unfavorable  effect  on  the  laminar  condition,  It 
assured  a  very  stable  free  stream  temperature  during  each  run.  This 
last  point  Is  very  Important  in  our  case  where  temperature  measurements 
were  also  taken. 

3.2  Flow  Measurements 

All  the  Information  collected  in  this  report  was  obtained  from 
three  basic  types  of  measurements:  wall  pressure,  pitot,  and  steady 
state  hot  wire  measurements.  The  following  sections  describe  all  three 
procedures  used  In  reducing  these  measurements. 

3.2.1  Wall  pressure  measurements,  location  of  separation  and 
reattachment  points 

The  wall  pressure  Is  an  important  parameter  in  this  type 
of  flow.  It  was  measured  as  stated  earlier  with  a  Decker  pressure  trans¬ 
ducer  and  the  result  reconverted  In  p.s.1.  using  the  transducer's 
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calibration  curve.  P  was  then  used  mostly  In  the  form  P/P^  where  P„ 

oo 

was  the  undisturbed  free  stream  condition.  The  location  of  separation 

and  reattachment  was  found  with  the  help  of  a  technique  used  by  Roshko 
14  1516 

and  l.iomke  and  others.  ’  It  consists  of  perturbing  the  flow  at  each 
pressure  tap  and  observing  whether  the  perturbation  gives  an  Increase  or 
a  decrease  In  the  measured  pressure.  If  the  perturbation  Is  downstream 
of  the  pressure  tap,  it  will  create  a  positive  pressure  gradient  and  a 
higher  pressure  will  be  recorded.  If  a  lower  pressure  Is  recorded,  It 
means  that  the  perturbation  Is  upstream  from  the  pressure  tap.  Experi¬ 
mentally  tire  state  was  realized  by  bringing  a  fine  wire  (diameter  of 
the  wire  *  diameter  of  the  pressure  tap  hole),  attached  to  the  arm  of  the 
traversing  mechanism,  to  one  side  of  each  pressure  tap  and  noting  the 
pressure  each  time.  The  same  procedure  was  then  repeated,  placing  the 
wire  at  the  other  side  of  the  pressure  taps,  and  finally  an  unperturbed 
pressure  was  recorded.  The  three  curves  were  plotted  on  the  same  graph 
and  the  point  where  they  Intersect  is  the  reattachment  point.  (At  that 
point,  where  the  flow  impinges  on  the  surface,  the  perturbed  pressure  and 
the  unperturbed  pressure  are  equal). 

This  technique  was  mainly  used  to  pick  up  the  reattach¬ 
ment  point,  the  separation  point  being  located  by  the  thin  oil  film 
technique.  (Using  the  fact  that  the  shear  after  separation  Is  negative 
and  ahead  of  separation  it  Is  positive,  the  oil  will  be  collected  by 
this  action  alcng  a  line  which  determines  the  separation  point.) 

Finally,  the  wall  pressure  measurements  were  used  to 
Investigate  the  effects  of  change  In  conditions  at  the  end  of  the  ramp. 
This  will  be  explained  more  extensively  in  the  next  chapter. 
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3.2.2  Pitot  measurements 

The  pitot  probe  described  earlier  was  used  to  take 
boundary  layer  stagnation  pressure  profiles  throughout  the  interaction 
region.  These  measurements  were  taken  for  two  reasons:  (a)  to  test 
the  hot  wire  data  procedure  described  in  the  next  section;  and  (b)  to 
use  Pt  (instead  of  wall  pressure)  along  with  the  hot  wire  data  to 
Investigate  the  presence  of  normal  pressure  gradients  in  the  boundary 
layer.  This  latter  quantity  is  difficult  to  measure  directly  due  to 
the  extremely  small  size  of  the  region  of  interaction. 

Figure  8  shows  a  typical  Mach  number  pitot  trace  ahead 
of  separation;  the  same  figure  also  shows  a  Mach  number  profile  computed 
from  hot  wire  data  and  wall  pressure.  Several  attempts  to  correct  the 
pitot  data  near  the  wall,  using  Homann's  approach,17  produced  the  same 
results  as  the  hot  wire.  Hot  wire  data  will  replace  the  pitot  results 
for  points  near  the  wall  in  all  that  follows.  This  means  that  when  the 
stagnation  pressure  of  the  pitot  is  used  with  the  hot  wire  measurements 
to  check  for  normal  pressure  gradients,  we  will  be  unable  to  detect  any 
pressure  gradients  in  the  bottom  part  of  the  boundary  layer  (roughly 
where  M  <  0.3).  It  is,  however,  very  unlikely  that  any  normal  pressure 
gradient  exists  tiiere,  because  as  we  approach  the  wall,  the  streamlines 
tend  to  become  parallel  to  it. 

3.2.3  Hot  wire  measurements 

Due  to  the  size  of  a  hot  wire  probe,  it  seemed  to  us 
the  most  practical  measuring  instrument  for  separated  flows.  Thus  it 
was  decided  that  this  would  be  the  main  source  of  data  for  our  investi¬ 
gation. 

The  hot  wire  measures  essentially  two  quantities:  the 
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resistance  of  equilibrium  (when  the  wire  is  not  heated)  and  the  heat 

loss  of  the  wire  when  It  Is  overheated  with  respect  to  the  flow  around 

it.  These  two  quantities  can  be  expressed  in  another  form;  namely,  the 

temperature  of  equilibrium  Te„  and  the  Nusselt  number  Nu.  Te  Is 

m  mm 

a  unique  function  of  rgm  and  the  electrical  resistivity  of  the  wire, 
and 


wSere  qm  is  the  heat  dissipated  in  the  wire  per  unit  area,  d  Its 
diameter,  kQ  the  thermal  conductivity  of  the  surroundings  at  the  stag¬ 
nation  conditions,  and  Tw^  the  temperature  of  the  heated  wire.  The 
subscript  m  stands  for  "measured." 

These  data  have  to  be  corrected  for  end  losses,  since  the  wire 
has  a  finite  length  and  is  soft  soldered  on  each  end  to  the  two  needle 
supports.  The  aspect  ratio  of  our  wire  was  of  the  order  of  300,  but,  as 

will  be  seen  later,  the  correction  for  end  losses  can  still  be  appreciable. 

18  19 

A  derivation  of  this  correction  is  given  by  Kovasny  and  Dewey,  and 
depends  essentially  on  a  parameter  S  given  by 


where  k  is  the  thermal  conductivity  of  the  wire  material  and  r  and 
rg  are  the  resistances  of  the  heated  and  adiabatic  wire.  This  quantity 
determines  two  correction  factors  CN  and  CR  which  give  us  the  Nusselt 
number  and  recovery  ratio  for  infinite  wires: 


and 
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"  '  CR  \ 

where  n  ■  yit  ;  %  =  Te^ 

Is  a  function  of  S  alone  and  is  a  function  of  S  and  the 
recovery  temperature  of  the  needle  supports.  Details  on  these  quanti¬ 
ties  are  given  in  the  Appendix. 

Now  that  our  measurements  are  corrected,  we  can  use 
the  available  Information  on  heat  loss  and  recovery  temperature  of 
infinite  cylinders  normal  to  the  stream.  Several  investigations  have 
shown  that  for  small  M  and  ReQ  (ReQ  =  pud/ut)  the  Nusselt  number 
is  a  function  of  these  two  quantities.  Dewey  has  formula,.ed  an  empirical 
equation  which  represents  a  curve  fit  of  the  existing  experimental  data. 
This  curve  fit  is  represented  by  an  equation  of  the  form: 

Nu0(RVM)  =  Nuo^Reo,ao)  ®(RVM) 

where  Nu0(ReQ,®)  represents  the  dependence  of  Nu^  on  ReQ  when 
M  >>  1  and  4>(ReQ,M)  is  the  departure  from  this  rriation  when  M  is 
no  longer  very  large.  These  relations  are  given  in  the  Appendix. 

If  we  call  n  the  recovery  temperature  T  /T^  we  know 
that  In  continuous  flow  n  is  a  unique  function  of  M;  it  varies  be¬ 
tween  1  for  M  =  0  and  0.95  for  M  >  2.  In  our  case,  however, 
using  very  low  free  stream  stagnation  pressures,  the  fluid  can  no  longer 
be  considered  a  continuum  at  the  scale  of  the  wire.  Dewey  has  proposed 
that,  if  we  write 


where  nc  is  the  value  of  n  at  the  continuum  limit  and  at  the 
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free  molecular  flow  limit,  then  n*  is  a  unique  function  of  the  Knudsen 

number  Kn^  for  a  given  M.  This  was  proved  by  the  data  of  several 

workers  in  the  field,  namely  Sherman,  Laufer  and  McClellan,  and  Dewey. 

The  relations  n  and  n?  as  a  function  of  M  and  n*  as  a  function 

of  Kn^  are  giver  in  the  Appendix. 

To  reduce  the  hot  wire  data  practically  requires  an 

iteration  scheme  between  the  end  loss  correction  and  the  heat  loss  and 

recovery  factor  relations.  This  computational  scheme  is  shown  later  in 

this  section.  First  we  will  explain  the  experimental  procedures. 

Across  the  two  needles  of  the  wire  support  was  mounted 

a  0.0002"  diameter  Wollaston  wire  ( Pt ,1 0%  Rh).  No  attempt  was  made  to 

measure  the  wire  diameter  directly  and  the  manufacturer's  specifications 

were  adopted.  The  wire  was  then  annealed  to  a  dull  glow  for  several 

minutes  and  brought  back  to  room  temperature  before  it  was  calibrated. 

This  procedure  was  carried  out  for  all  wires  to  avoid  a  great  change  in 

the  electrical  resistivity  coefficient  after  each  heating.  The  wire  was 

then  calibrated  and  the  coefficients  a  and  r  of  a  relation  of  the 

r 

form  T  =  ar  +  rr  was  computed.  (In  our  case  the  range  of  temperature 
is  small  and  a  is  a  constant  independent  of  T).  The  resistance  per 
foot  at  20°C  was  given  by  the  manufacturer.  Using  this  value,  the  length 
of  the  wire  i  was  computed  (£w  was  of  the  order  of  5.4  x  10  ft). 

The  wire  was  then  placed  in  the  test  section,  and  at  each  point  of  the 
boundary  layer,  two  values  of  overheat  coefficients  [overheat  =  (r  -r  )/r  ] 

W  t?  I? 

smaller  than  0.08  were  taken  and  the  voltage  drop  across  the  wire  was  re¬ 
corded  (see  previous  section).  The  resistance  was  then  plotted  against 
i  and  the  slope  computed  (actually  this  was  calculated  directly  without 
plotting  because  1  is  a  linear  function  of  *■  for  small  overheats). 


With  the  calibration  curve  of  the  wire  and  r  ,  Te„  was  then  computed. 

em  m  r 

The  following  computational  scheme  was  programmed  and 
computed  on  the  CDC  6400  of  the  Computer  Center  at  the  University  of 
California  at  Berkeley. 

1.  Assume  =  1  (i.e.,  Tem  =  Tt).  With  this  value  of 
Te  compute  kQ  and  Num.  Compute  S  and  hence  NuQ,  and  from  the 
relation  Nu0(ReQ,M)  find  M  [note  that  Nuo(ReQ,M)  can  also  be 
written  Nuq(M,  T^,  P) ] . 

2.  Once  M  is  known,  and  using  the  same  value  for  T^, 

compute  Kn^  with  either  (pitot)  or  P  (wall  pressure).  Also 

compute  and  nc.  and  hence  find  n* 

3.  With  the  value  found  for  M  compute  ns  (ns  =  Ts/Tq 
support  recovery  ratio),  and  hence  find  the  end  loss  recovery  ratio 
CR  (see  Appendix). 

4.  Knowing  CR  and  rj  compute  r^.  Go  to  1  and  use 

the  computed  value  for  nm  instead  of  =  1. 

This  iteration  converges  very  rapidly  and  usually  no 
changes  to  the  fourth  decimal  place  were  observed  after  the  third  iter¬ 
ation.  In  the  above  computation,  a  Sutherland's  viscosity  law  was  used. 

The  knowledge  of  M,  T^,  and  either  the  wall  pressure 
P  or  the  pitot  pressure  P^  enabled  us  to  compute  all  thermodynamic 
and  dynamic  quantities  in  the  boundary  layer. 
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4.0  RESULTS  AND  DISCUSSION 


The  data,  collected  and  reduced  by  the  methods  described  In  the  pre¬ 
vious  sections,  will  now  be  analyzed. 

We  will  first  give  a  general  description  of  these  results  and  attempt 
to  relate  them  to  previous  experimental  and  theoretical  papers  treating 
this  problem.  In  the  second  section  we  will  present  a  detailed  analysis 
of  the  data  and  attempt  to  explain  some  aspects  of  this  complicated  flow 
in  the  light  of  the  measurements  of  local  boundary  layer  quantities  which 
were  obtained. 

4.1  General  Discussion.  Comparison  with  Previous  Results 

The  first  measurements  obtained  were  the  pressure  distributions 
in  the  streamwise  direction  for  different  free  stream  conditions  and  com¬ 
pression  angles  a.  These  measurements  show  very  clearly  a  pressure  rise 
ahead  of  the  corner,  then  a  constant  pressure  region  that  extends  to  a 
small  distance  downstream  of  the  corner,  followed  by  another  pressure  rise 
to  reattachment  and  beyond.  (See  Fig.  5.)  This  compares  qualitatively 
very  well  with  previous  results  of  Chapman  et  al.,  Lewis,  etc.  Measuring 
the  pressure  at  separation  p<.  and  the  pressure  of  the  plateau  region  p 

O  p 

(the  separation  point  being  located  by  the  method  described  earlier),  we 
plotted  In  Fig.  9 


Pp  -  Pp 
Pp  ■'w 


and 


Ps  -  Pp 
Pp  -'W 


versus  M.  This  correlation  in  which  =  c^ 
suggested  earlier  (see  Ref.  1  for  more  details 


/[(cf  )  J]  has  been 
o  To  R*10° 

).  On  the  same  figure  we 


also  used  results  from  Chapman's  Investigations  (open  symbols)  which 
proved  the  quantitative  agreement  of  our  results  with  his,  at  least  up 
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to  the  plateau  pressure.  The  pressure  distributions  i’i  the  second  recom¬ 
pression  region  could  not  be  compared  qualitatively  with  any  result  because 
there  is  no  way  to  correlate  the  measurements  in  that  region;  more  details 
concerning  this  point  will  be  pointed  out  in  the  next  section. 

The  locations  of  the  points  of  separation  and  reattachment  also 
proved  to  be  comparable  to  those  found  by  Chapman.  The  pressure  at  the 
separation  point  was  always  slightly  mere  than  half  of  the  plateau 
pressure.  The  pressure  at  reattachment  cannot  be  given  in  terms  of  the 
difference  between  the  final  pressure  and  the  plateau  pressure;  in  fact, 
it  is  wrong  to  give  such  a  value,  for,  as  shown  in  11,  and  as  will  also  be 
shown  in  the  next  section,  the  final  pressure  downstream  of  the  reattach¬ 
ment  point  can  have  different  values  without  affecting  the  main  flow. 

The  hot  wire  data  complemented  by  either  the  wall  pressure  or 
the  pitot  total  pressure  can  give  us  the  values  of  all  the  thermodynamic 
and  dynamic  quantities  in  the  boundary  layer.  Using  the  data  reduction 
procedure  explained  in  Section  3.2.3,  profiles  similar  to  those  in  Figs. 

10  to  13  were  obtained.  A  total  of  20  such  profiles  were  obtained  at 
different  streamwise  locations  for  two  compression  angles  and  free  stream 
conditions.  It  is  difficult  to  compare  these  results  with  any  previous 
work,  for,  to  our  knowledge,  no  measurements  of  this  kind  have  been  per¬ 
formed.  The  Mach  number  profiles  across  the  boundary  layer  are  quali¬ 
tatively  comparable  to  those  obtained  by  a  pitot  in  (10);  however,  our 
measurements  seem  to  be  more  accurate  in  the  region  close  to  the  wall. 
Figure  10  shows  velocity  profiles  at  two  stations  ahead  of  separations. 
Their  deformation  from  a  Blasius  type  profile  to  a  near  separation 
profile  is  quite  clear.  Figure  11  shows  the  static  and  stagnation  en¬ 
thalpy  profiles  corresponding  to  the  previous  figure.  These  results 
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seem  to  agree  qualitatively  very  well  with  theoretical  predictions.  Note 
that  the  stagnation  enthalpy  is  not  constant  across  the  boundary  layer, 
but  goes  from  a  value  inferior  to  the  free  stream  stagnation  enthalpy, 
increases  with  y,  overshoots  the  free  stream  value  and  decreases  back 
to  this  value  at  the  edge  of  the  boundary  layer.  Note  also  that  the 
thermal  boundary  layer  is  thicker  than  the  velocity  boundary  layer,  al¬ 
though  by  a  small  amount.  This  obviously  comes  from  the  fact  that  the 
Prandtl  number  is  smaller  than  one. 

Similar  remarks  can  be  made  on  the  profiles  of  Figs.  11  and  12; 
we  can  add,  however,  in  the  case  of  Fig.  12  that  there  is  clear  evidence 
there  of  the  existence  of  a  reversed  flow.  In  this  figure  the  minimum 
value  of  u/u'  corresponds  to  the  zero  velocity  line;  the  values  of 
u/u 1  below  this  point  correspond  to  the  reversed  flow,  but  the  hot  wire 
response  is  independent  of  the  stream  direction. 

We  give  these  results  in  this  section  to  draw  a  general  picture 
of  the  flow  confirming  the  physical  qualitative  “"masoning  proposed  In  the 
introduction.  A  thorough  analysis  of  these  daf  is  carried  out  in  the 
next  section. 

4.2  Detailed  Study 

4.2.1  Existence  of  normal  pressure  gradients 

Although  the  assumption  of  zero  normal  pressure  gradient 
(dp/dy  =  0)  is  always  used  by  theoreticians  attempting  to  solve  this 
problem,  no  final  experimental  proof  has  been  advanced  which  either  proves 
or  disproves  this  basic  hypothesis.  As  we  mentioned  in  earlier  sections, 
our  experiments  included  the  indirect  computations  of  p  across  the 
boundary  layer;  the  results  are  shown  in  Figs.  14a  and  14b.  Before  dis¬ 
cussing  these  results,  let  us  make  some  pertinent  remarks  to  avoid 
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misinterpretations  due  to  large  errors  caused  by  experimental  uncertain¬ 
ties.  To  compute  p  the  following  parameters  were  measured: 

Re  .  Num  with  hot  wire  probe 
m  m  r 

Pt  with  the  pitot  probe. 

However,  it  was  impossible  to  measure  them  simultaneously  (the  traversing 
mechanism  can  hold  only  one  probe).  So  they  were  measured  on  two  differeit 
runs,  taking  the  following  precautions: 

--  The  free  stream  conditions  are  identical  in  the  two  runs. 

This  was  relatively  easy  to  accomplish  since  the  tunnel 
was  mostly  used  at  its  lowest  stagnation  pressures. 

—  The  position  of  the  probes  with  respect  to  the  wall  was 
carefully  measured  using  magnifying  optics. 

The  first  precaution  proved  very  effective  because  by  taking  wall  pressure 
measurements  on  the  two  runs  they  fall  within  2.5%  from  each  other.  In 
other  words,  the  flow  being  measured  by  the  two  instruments  at  different 
runs  was  the  same.  The  second  precaution  was  not  fully  satisfactory;  if 
an  error  of  0.002"  is  committed  on  the  distance  between  the  probe  tip  and 
the  wall  this  can  induce  errors  on  p^  of  the  order  of  30%  in  the  region 
of  large  shear.  This  is  due  to  the  fact  that  in  this  region  all  our 
measured  quantities  (Rem,  Num,  P^.)  vary  rapidly  in  the  y  direction. 

Note  that  for  this  very  same  reason,  the  error  will  be  smaller  in  the 
regions  of  low  shear,  namely  near  the  wall  (separated  profiles),  and 
around  the  edge  of  the  boundary  layer.  Also,  we  will  remember  (see  Section 
3.2.2)  that  the  pitot  data  near  the  wall  (M  <  0.3)  was  corrected  by  using 
the  hot  wire  data  along  with  wall  pressure.  This  means  that  no  pressure 
gradients  are  measurable  there  (using  our  method). 
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Figures  14a ,b  show  plots  of  P/Pwa^  as  functions  of 
y  at  different  streamwlse  locations.  A  first  remark  pertaining  to  what 
was  said  above  is  that  P/Pwa^  can  be  both  *  than  1  Inside  the  bound¬ 
ary  layer.  This  is  physically  difficult  to  Interpret  and  may  very  well 
be  due  solely  to  error  contribution. 

These  plots  show  a  clear  tendency  toward  negative  (very 
small)  normal  pressure  gradients  near  the  edge  of  the  boundary  layer  mainly 
at  the  stations  where  dp/dx  is  Urge  (near  separation  and  reattachment, 
but  not  in  the  plateau  region).  TMs  could  be  interpreted  by  the  fact 
that,  moving  the  probes  normal  to  tne  wall,  we  will  be  crossing  different 
compression  lines  (in  the  compression  waves)  that  are  present  In  the  super¬ 
sonic  upper  layer  of  the  boundary  layer.  However,  it  can  be  safely  said 
that  no  substantial  normal  pressure  gradients  are  present  In  the  flow  we 
investigated  (compression  angle  a  =  11°). 

Another  tist.  based  on  the  previous  data  was  to  plot  the 
wall  pressure  and  the  pressure  at  the  edge  of  the  boundary  layer  versus 
x/JL  Figure  15  shows  such  a  plot  and  confirms  the  conclusions  drawn  above. 

On  this  same  figure  we  also  shnw  the  pressure  distribu¬ 
tion  as  computed  f>om  the  displacement  thickness  profile,  using  Prandtl 
isentropis  compression  relationships  with  the  turning  angle  A0  ■  d62/dx. 
This  computed  free  stream  pressure  agrees  very  well  with  the  previous 
values. 

It  is  worth  noting  that  the  wall  pressure  downstream  of 
reattachment  always  reaches  values  larger  than  the  invlscld  pressure  (down¬ 
stream  of  an  oblique  shock  originating  at  the  corner).  This  is  not  true 
for  the  computed  free  stream  pressures.  Although  the  difference  between 
the  two  values  is  very  small,  it  is  quite  probable  that  if  the  compression 


angle  is  made  larger,  non-negl igible  normal  pressure  gradients  mir,ht  exist 
there.  For  this  reason  the  only  conclusion  we  can  safely  propose  is  that, 
provided  the  compression  angle  is  not  very  large  and  the  boundary  layer  is 
laminar,  no  evidence  of  the  existence  of  large  normal  pressure  gradients 
is  observed. 

4.2.2  Static  total  enthalpy;  wall  recovery  ratio 

Variations  of  static  enthalpy  in  the  y  direction  are 
represented  In  Figs.  16  and  17.  Ahead  of  the  interaction  h  profiles 
agree  very  well  qualitatively  with  known  data  for  a  flat  plate.  As  the 
interaction  starts,  a  nearly  constant  temperature  layer  starts  to  develop 
near  the  wall.  Moving  further  downstream  this  layer  becomes  still  thicker, 
varying  roughly  like  the  thickness  of  the  separated  bubble.  The  enthalpy 
distribution  in  the  region  of  high  shear  is  still  similar  to  that  of  an 
attached  boundary  layer. 

As  the  flow  approaches  reattachment,  the  inverse  trend 
is  observed;  namely,  the  region  of  nearly  constant  enthalpy  becomes  thinner 
while  the  variation  of  enthalpy  in  the  high  shear  region  is  still  qualita¬ 
tively  similar  to  that  of  an  attached  laminar  boundary  layer.  This  trend 
continues  after  reattachment,  but  the  changes  become  smaller  as  we  move 
downstream  for  the  case  of  a  laminar  boundary  layer.  If  the  boundary  layer 
begins  to  become  turbulent,  the  enthalpy  profile  becomes  fuller.  This,  by 
the  way,  is  also  observed  in  the  shadowgraph  when  the  bright  region  around 
the  edge  of  the  boundary  layer  begins  to  disappear.  This  effect  was 
explained  in  a  previous  section. 

A  closer  look  at  the  separated  bubble  reveals  that  the 
enthalpy  does  not  vary  in  the  y  direction  up  to  about  halfway  between 
the  zero  velocity  line  and  the  separating  streamline.  Furthermore,  the 
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enthalpy  does  not  vary  appreciably  In  the  x  direction  In  a  region  extend¬ 
ing  from  the  beginning  to  the  end  of  the  plateau  constant  pressure  region. 
(Note  that  the  separated  bubble  extends  both  upstream  and  downstream  of  the 
constant  pressure  plateau  region). 

The  previous  remarks  suggested  a  comparison  between  the 
adiabatic  wall  temperature  and  the  equilibrium  temperature  of  a  flat  plate 
in  a  compressible  flow  at  the  same  local  free  stream  Mach  number.  Using 
the  value  0.723  for  Prandtl  number,  the  theoretical  prediction  (see  Ref.  21) 
was  computed  at  different  stations  downstream  of  the  leading  edge,  using  at 
each  station  the  free  stream  Mach  number  measured  with  the  hot  wire.  Ex¬ 
trapolating  the  temperature  profile  to  the  wall,  we  find  j  /T •  at  differ- 

w 

ent  stations.  Figures  :C  and  19  show  two  plots  of  T  /T '  compared  with 
the  theoretical  predictions  of  Ref.  21. 

The  agreement  is  good  in  general.  There  is,  however,  a 
consistent  tendency  for  a  higher  measured  wall  temperature  even  before 
the  interaction  starts.  The  error  may  be  due  largely  to  the  conduction  of 
heat  inside  the  metallic  model  from  regions  of  high  temperatures  to  cooler 
regions . 

It  is  safe,  however,  to  assume  that  the  temperature  inside 
the  separated  bubble  is  equal  to  the  temperature  that  the  wall  would  take 
if  the  boundary  layer  remained  attached  and  the  local  free  stream  Mach 
number  is  that  of  the  separated  flow. 

4.2.3  Shear  stress 

The  shear  stress  y(3u/3y)  was  calculated  using  Sutherland's 
viscosity  law  for  y  and  graphically  determining  3u/3y  from  the  velocity 
profiles  mentioned  earlier.  A  graphical  interpolation  was  used  when  com¬ 
puting  the  shearing  stress  on  the  zero  velocity  line.  This  procedure  is 


based  on  the  fact  that  u  is  smooth  on  the  zero  velocity  line  and  au/9y 
approximated  by  the  tangent  has  the  same  value  on  each  side  of  this  point. 
Due  to  the  inaccuracy  in  determining  u  for  the  low  velocities  in  the 
separated  bubble,  the  procedure  described  inevitably  induces  large  errors. 
The  following  results  should  therefore  be  regarded  as  more  descriptive 
than  quantitative. 

Figures  20  and  20  show  the  variation  of 
[u( 9u/3y)]r 

with  x.  The  reference  value  [u(3u/3y)]r  is  taken  to  be  the  shear  stress 
ahead  of  interaction  (x/l  =  0.25).  At  this  station  the  boundary  layer 
profile  is  essentially  a  Blasius  profile. 

The  variation  of  this  quantity  shows  the  fast  decrease  in 
the  wall  shear  stress  ahead  of  separation  and  after  reattachment.  In  the 
separated  bubble,  the  changes  are  much  smaller  and  a  clear  maximum  is 
observed  near  the  corner. 

The  increase  of  the  shear  after  reattachment  is  different 
in  t.iese  two  figures.  Figure  20  shows  a  pure  laminar  reattachment  (where 
the  boundary  layer  stays  laminar  several  boundary  layer  thicknesses  down¬ 
stream  of  reattachment),  ai.d  in  Fig.  21  the  flow  becomes  transitional 
slightly  downstream  of  reattachment.  This  part  of  the  flow  will  be  dis¬ 
cussed  in  more  detail  in  the  last  section. 

It  is  clear  from  the  velocity  profiles  in  the  separated 
bubble  that  the  shear  stress  at  the  wall  is  almost  equal  but  opposite  to 
the  shear  on  the  zero  velocity  line.  (We  used  the  latter  in  our  graphs 
because  the  error  in  computing  it  seemed  smaller  than  that  which  could 
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result  from  the  uncertainty  in  the  relative  positions  of  the  probe  with 
respect  to  the  wall.)  Consequently,  we  can  make  a  qualitative  comparison 
of  our  results  with  those  of  Holt,  which  seems  to  confirm  the  existence 
of  a  minimum  (negative  maximum)  of  the  shear  stress  near  the  corner. 

4.2.4  Flow  Inside  the  separated  bubble 

At  different  locations  in  the  streamwise  direction,  the 

zero  velocity  (u  =  0)  point  was  determined  as  described  in  the  previous 

section.  This  allows  us  to  find  at  the  same  x  the  point  A  on  the  mass 

flux  profile  which  belongs  to  the  separating  streamline 

u=0  point  A 

/  pu  dy  =  /  pu  dy 

o  u=o  point 

This  simply  says  that  the  mass  flux  between  the  u  =  0  line  and  the  wall 
is  equal  to  the  mass  flux  between  the  u  -  0  line  and  the  separating 
streamline  at  the  same  streamwise  position.  Figure  22  illustrates  the 
above  integration  procedure  which  was  carried  out  graphically  by  an  Amsler 
integrater.  A  similar  procedure  was  used  to  graph  streamlines  inside  the 
separated  bubble. 

Figure  23  shows  the  edge  of  the  boundary  layer,  the  dis¬ 
placement  thickness  and  several  streamlines  (numbered  from  1  to  5)  inside 
the  separated  bubble. 

Defining 

—  1  *5 
pu  (y  =  0,  y  =  y,)  =  —  /  pu  dy 

5  y5  o 
and 

—  1  Z 

pv  (x  =  x  ,  x  =  l)  =  t—  /  pv  dx 

5  s  x$ 

The  following  values  give  an  idea  of  the  magnitude  of 
the  quantities  of  Interest  inside  the  separated  bubble: 
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at  x/ l  »  1 


pu  (y  =  0,  y  =  y,) 
- — “ 


0.028 


pu  (y  =  y5»  separating  streamline) 


p'u1 


0.067 


and  along  the  u  =  0  line 

pv  (x  ■  x  ,  x  ■  t) 

- TT-i -  =  0.0018 

P  U 

pv  (x  =  l,  X  =  Xp) 

- -T-, - —  =  0.0022 

P  u 

As  for  the  velocities,  it  is  clear  from  Fig.  23  that  the 
largest  velocities  in  the  separated  bubble  are  along  the  separating  stream¬ 
line.  In  the  reversed  flow  itself,  the  order  of  magnitude  of  u  was  3% 
to  10%  of  the  free  stream  velocity  and  that  of  v  (along  the  u  =  0  line) 
above  0.3%  to  0.7%  of  u‘ . 

Figure  23  also  includes  points  near  the  separating  stream¬ 
line  that  were  computed  such  that  the  total  pressure  on  these  points  is 
equal  to  the  wall  pressure  at  the  reattachment  point.  The  fact  that  these 
points  are  quite  near  the  separating  streamline  means  that  Chapman's 
isentropic  recompression  hypothesis  is  a  valid  assumption.  A  double  check 
concerning  this  point  is  also  possible.  Along  the  separating  streamline 
we  have: 

p 

d  ,  u  \  _  dp  ,  9t 

3*  (  “T  >  '  *  Hx  97 


Our  result  shows  that  3x/3y  is  small  on  the  separating  streamline  (the 
curvature  of  the  velocity  profile  is  very  large  there).  9T/3y  increases 
when  we  get  near  the  reattachment  but  so  does  dp/dx,  and  it  is  not 
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surprising  that  the  isentropic  recompression  hypothesis  is  a  very  good 
approximation. 

It  is  worth  noting  here  that  the  accuracy  of  the  hot  wire 
measurements  decreases  substantially  when  the  velocity  becomes  small.  The 
above  results  should  then  be  considered  as  descriptibe,  giving  only  an 
order  of  magnitude  of  the  quantities  of  interest.  Again  it  must  be  pointed 
out  that  the  verification  of  Chapman's  hypothesis  is  valid  provided  the 
angle  of  compression  a  is  small  (a  <  13°),  and  that  the  boundary  layer 
is  laminar. 

4.2.5  Detailed  study  of  the  reattachment  region 

Relatively  little  is  known  about  the  mechanism  of  reattach¬ 
ment  as  compared  with  the  separation  process  where  the  free  interaction 
theory  seems  so  successful.  It  seems  to  us  that  the  main  experimental 
difficulty  encountered  in  that  region  is  the  fact  that  the  boundary  layer 
often  becomes  turbulent  when  it  experiences  the  negative  pressure  gradient 
near  recompression.  In  the  following  study  a  special  emphasis  is  put  on 
this  point  for,  as  we  will  see,  a  transitional  reattachment  can  be  radi¬ 
cally  different  from  a  laminar  one. 

As  stated  earlier.  Model  B,  which  was  used  for  this  study, 
nad  a  provision  for  changing  the  conditions  downstream  of  the  reattachment 
print.  More  precisely,  we  could,  at  a  chosen  point  on  the  ramp,  create  a 
positive  or  negative  pressure  gradient  of  any  intensity.  Furthermore,  this 
pressure  gradient  could  be  imposed  either  suddenly  by  a  sharp  deformation 
of  tne  wall  (deformation  spread  over  a  distance  of  the  order  of  one  boundary 
layer  thickness),  or  gradually  (deforming  the  wall  over  a  distance  of 
several  boundary  layer  thicknesses). 


Performing  the  experiment  described  above  at  different 
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locations  downstream  of  R,  we  found  as  in  Ref.  11  that  there  exists  a 

point  C,  downstream  of  which  the  conditions  of  the  flow  do  not  affect 

the  separated  flow  and  the  reattachment  upstream  of  C^.  The  pressure 

distribution  downstream  of  C,  can  lie  between  two  limits  P  and 

1  max 


Pm<n;  the  range  P 


max  min 


P  .  increases  moving  downstream  (see  Fig.  24). 


Disturbing  the  flow  ahead  of  C.  usually  changes  by 
more  than  5%  the  pressure  distribution  ahead  of  reattachment.  However, 
a  gradual  pressure  gradient  could  be  applied  at  some  location  without 
disturbing  the  flow  upstream  of  this  location.  By  a  gradual  pressure 
gradient  we  mean  a  generally  small  pressure  change  spread  over  a  distance 
larger  than  three  or  four  boundary  layer  thicknesses.  This  can  be  achieved 
up  to  a  point  slightly  downstream  of  R.  We  call  this  point  C^,  the 
second  critical  point.  Upstream  of  any  change  of  the  pressure  dis¬ 
tribution  however  small  or  gradual  changes  the  pressure  distribution 
upstream  by  more  than  5%.  This  result  has  been  observed  earlier,  and  we 
seem  to  have  the  same  problem  locating  the  second  critical  point  precisely. 
This  may  be  largely  due  to  the  fact  that  near  the  reattachment  point  the 
velocities  near  the  wall  are  extremely  small. 

It  was  noted  that  C£  is  usually  located  at  a  point  on 
2  2  .  . 

the  pressure  distribution  where  d  p/dx  starts  to  become  negative;  i.e., 
near  the  end  of  the  almost  linear  increase  of  pressure  with  x  (see  Figs. 
24  and  25).  As  in  Ref.  11  the  point  was  found  to  have  an  interesting 
property.  Writing 


p  the  pressure  at  c? 

c2 

Pp  the  plateau  pressure 

x„  the  location  of  c9 

c2  2 
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x  the  location  of  the  end  of  the  plateau  constant 
p  pressure  region 

and  plotting  (p  -p)/(p  -p  )  versus  (x  -x)/(x  -x  )  (Fig.  26),  we 

P  c2  l2  ep 

see  that  we  can  correlate  the  pressure  distributions  found  at  different 
free  stream  conditions  up  to  C£.  This  figure  also  includes  data  from 
previous  investigations  where  we  estimated  the  location  of  using 
the  property  mentioned  previously. 

Noting  that  the  pressure  distribution  slightly  downstream  of 
the  plateau  pressure  increases  linearly  up  to  C^,  we  ask  ourselves 
if  and  how  this  slope  could  be  related  to  the  free  stream  conditions  and 
model  geometry.  In  order  to  compare  data  from  different  sources  Instead 
of  look.,  jt  p  versus  x/£  distribution,  it  is  more  meaningful  to 
look  at  p  versus  Re^.  The  most  obvious  parameter  on  which  the  slope 
of  p  might  depend  is  the  angle  at  which  the  separated  free  shear  layer 
impinges  on  the  ramp.  Call  this  angle  (a-B).  where  B  Is  the  angle 
of  the  free  shear  layer  with  respect  to  the  flat  plate. 

Figure  32  shows  the  dependence  of  the  slope  of  p  versus  Rex 
on  (a-8)  (some  of  the  points  were  taken  from  other  references  for  larger 
a  and  M).  This  figure  suggests  a  rather  simple  dependence;  however,  it 
should  not  be  regarded  as  an  accurate  correlation  for  the  pressure  distri¬ 
bution  near  the  reattachment  point.  It  should  be  noted  here  that  B  Is 

a  function  of  the  plateau  pressure,  which  is  itself  a  function  of  M,Re  , 

o 

so  that  the  quantity  (a-B)  includes  both  the  free  stream  conditions  and 
model  geometry.  The  points  of  Fig.  32  that  are  computed  from  Ref.  1 
correspond  to  leading  edge  separation;  the  fact  that  they  fit  well  with 
the  other  points  suggests  that  the  mechanism  of  reattachment  does  not 
depend  strongly  on  the  mechanism  that  induced  separation,  but  rather  on  the 
mixing  process  after  separation. 


Although  the  location  of  could  not  be  given  with  an  accu¬ 

racy  of  better  than  plus  or  minus  one  boundary  layer  thickness,  it  was 
consistently  found  that  its  distance  from  the  reattachment  point  did  not 
change  very  much  with  small  changes  in  the  free  stream  Reynolds  number. 
Such  was  not  the  case  for  ;  it  was  observed  that  a  small  change  of 
ReQ  could  move  several  boundary  layer  thicknesses  even  when  the 
pressure  distribution  in  the  separated  flow  did  not  change  by  more  than 
3%.  This  casts  some  doubts  about  the  critical  importance  of  such  a 
point  on  the  reattachment  process.  Furthermore,  it  was  noticed  that 
moves  downstream  with  increasing  ReQ.  This  suggested  a  relation  of  this 
point  with  the  occurrence  of  transition. 

Figure  27  shows  velocity  profiles  downstream  of  the  reattach¬ 
ment  point.  Note  that  the  boundary  layer  does  not  change  appreciably 
over  a  distance  of  about  ten  boundary  layer  thicknesses.  In  the  flow  that 
corresponds  to  these  profiles,  the  critical  point  was  not  detectable 
before  the  end  of  the  ramp. 

Figure  28  shows  transitional  velocity  profiles  downstream  of 
the  reattachment  point.  The  flow  upstream  of  the  second  critical  point 
was  identical  to  that  corresponding  to  Fig.  26.  However,  in  this  new 
case  the  first  critical  point  was  located  and  had  the  properties  described 
earlier.  The  last  two  figures  suggest  that  the  physical  meaning  of  this 
point  is  that  the  boundary  layer  has  picked  up  enough  momentum  to  enable 
it  to  overcome  finite  sudden  pressure  gradients.  Figures  29  and  30 
further  support  this  argument.  Looking  at  the  relative  slopes  of  the 
displacement  thickness  with  boundary  layer  edge,  we  see  that  in  a  transi¬ 
tional  flow  there  is  substantially  more  momentum  entering  the  boundary 
layer.  The  location  of  the  first  critical  point  could  thus  be  related 
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to  the  occurrence  of  turbulence. 

Attempts  to  relate  the  location  of  the  first  critical  point  in 
terms  of  some  parameter  involving  mean  quantities  of  the  boundary  layer 
yielded  Fig.  31.  In  this  figure  we  chose  the  mean  Mach  number  defined  as: 

i  6 

H  «  r  /  M  dy 
4  0 

as  this  parameter,  and  we  plotted  M  versus  x/J l  for  laminar  and  transi¬ 
tional  reattachment.  Even  though  R  becomes  larger  than  one  in  the 
transitional  case,  it  is  clearly  not  the  parameter  to  distinguish  between 
the  subcritical  and  supercritical  conditions  of  the  Crocco  Lees,  Lees  and 
Reeves  theory.  The  theoretical  parameters  of  these  theories  would  indi¬ 
cate  a  subcritical  condition  for  an  adiabatic  Blasius  profile.  This  is 
not  the  case  of  M  ahead  of  separation  in  Fig.  31. 

The  fact  that  the  first  critical  point  was  not  detectable  before 
the  end  of  the  ramp  in  the  case  of  pure  laminar  reattachment,  does  not 
necessarily  imply  the  non-existence  of  this  point.  However,  if  this  point 
exists  it  will  be  at  several  score  boundary  layer  thicknesses  from  the 
reattachment.  Furthermore,  there  is  some  doubt  as  to  the  possibility  of 
a  laminar  reattaching  boundary  layer  becoming  "supercritical"  (in  the 
sense  of  the  theories  mentioned  above)  because  when  the  boundary  layer  is 
fully  rehabilitated,  it  will  tend  to  become  a  Blasius  profile,  and  this 
profile  is  known  to  be  subcritical.  This  remark  is  clearly  not  valid  for 
a  laminar  wake  flow  where  the  mean  properties  of  the  boundary  layer  after 
the  rear  stagnation  point  will  eventually  imply  a  supersonic  behavior. 

This  will  also  be  the  case  for  a  turbulent  flow  where  the  distance  over 
which  the  boundary  layer  receives  enough  momentum  to  rehabilitate  itself 
will  be  substantially  smaller  and  close  to  C2. 
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Figures  29  and  30  also  show  a  clear  neck  in  the  boundary  layer 
thickness  when  the  flow  is  transitional;  this  is  not  observed  for  pure 
laminar  reattachment  where  the  boundary  layer  decreases  very  slowly  after 
the  second  critical  point  and  reaches  an  almost  constant  thickness  for 
boundary  layer  thicknesses  downstream. 

To  state  briefly  the  results  of  this  section,  we  can  say: 

1.  There  exists  a  critical  point  C2  slightly  downstream  of  reattach¬ 
ment  which  plays  an  important  role  in  the  mechanism  of  reattachment  (see 
Fig.  26)  and  whick  marks  the  beginning  of  the  rehabilitation  of  the  re¬ 
attached  boundary  layer. 

2.  If  transition  to  turbulent  flow  occurs  downstream  of  C^,  another 
critical  point  C-|  can  be  located.  This  point  has  the  property  that  down¬ 
stream  of  it  the  conditions  do  not  affect  the  flow  upstream,  provided  they 
lie  in  some  range.  Cj  ca  mot  be  located  for  laminar  flow  (at  least  close 
to  R),  and  there  is  doubt  about  its  theoretical  importance  for  reattaching 
laminar  boundary  layers. 
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5.0  CONCLUSIONS  AND  SUGGESTIONS  FOR  FUTURE  WORK 

A  detailed  experimental  study  has  been  made  of  the  two--imensional 
laminar  flow  over  a  compression  corner  with  Mach  numbers  around  2.5.  The 
following  conclusions  were  reached. 

1.  For  a  small  compression  angle  a  {a  <  13°)  and  a  laminar  boundary 
layer  there  is  no  evidence  of  appreciable  normal  pressure  gradients  in  the 
boundary  layer  ahead,  during,  and  after  interaction. 

2.  The  separated  bubble  has  sensibly  constant  enthalpy  equivalent  to 
the  wall  enthalpy  of  a  flat  plate  with  attached  boundary  layer  and  an  equiva¬ 
lent  free  stream  Mach  number. 

3.  The  shear  stress  along  the  u  =  0  line  is  usually  not  larger  than 
10%  of  the  shear  stress  at  the  wall  ahead  of  the  interaction.  The  maximum 
of  the  shear  stress  along  the  u  =  0  line  occurs  near,  but  slightly  down¬ 
stream  of  the  corner. 

4.  Evidence  of  a  reversed  flow  about  10%  the  velocity  of  the  free 
stream  exists.  The  maximum  of  u  in  the  seoarated  bubble  is  along  the 
separating  streamline,  and  is  always  smaller  than  half  the  local  speed  of 
sound.  Provided  a  is  small  and  the  boundary  layer  is  laminar,  the 
assumption  of  isentropic  recompression  is  a  good  approximation. 

5.  Two  critical  points  were  located  J  •*  the  separated  region.  One 
of  them  seems  to  be  greatly  affected  by  the  occurrence  of  transition  and 
is  not  detectable  for  pure  laminar  flow  over  a  distance  of  more  than  15 
boundary  layer  thicknesses  downstream  of  reattachment.  The  other  is  con¬ 
sistently  near  the  reattachment  point  and  seems  to  play  an  important  role 
in  the  reattachment  recompression  process. 

Due  to  the  limitations  mentioned  earlier,  some  aspects  of  the  problem 
could  not  be  investigated.  Future  work  based  on  the  same  techniques  could 
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clarify  the  following  points: 

1.  Existence  of  normal  pressure  gradients  and  validity  of  isentroplc 
recompression  along  dividing  streamline  when  a  is  no  longer  small  and/or 
when  the  boundary  layer  is  turbulent. 

2.  Attempts  to  locate  the  first  critical  point  with  a  longer  ramp. 
Test  the  critical  Importance  of  this  point,  if  any. 

3.  Detailed  enthalpy  study  for  the  case  of  highly  cooled  wall. 
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APPENDIX 


1. 


End  Losses  Correction  Factors 

The  Nusselt  number  end  loss  correction  factor  Is  given  by: 
_  * 


CN 


1  +  a 
1  ♦  a 


w 

w 


_  * 

where  a^/a  w  as  a  function  of  S  Is  represented  In  Fig.  A-l. 

The  recovery  temperature  correction  factor  Is  given  by 


CR  •  [1  -«<£)]  (1  -  u,)*1 

TTI 

where  uj/(  1  -cu)  as  a  function  of  S  Is  given  In  Fig.  A-l. 

A  series  of  runs  was  made  to  determine  the  needle  support 
temperature  T$.  A  resistance  thermocouple  (actually  using  a  hot  wire 
as  a  resistance  thermocouple)  was  mounted  on  one  needle  and  the  temper- 
at  different  Mach  numbers  was  measured.  The  results  are  shown  In  Fig. 
A-2,  together  with  the  predicted  recovery  temperature  of  cones  In  super¬ 
sonic  flow  using  laminar  boundary  layer  theory.  It  Is  this  predicted 
recovery  temperature  which  was  used  In  the  computation  of  n$.  Although 
the  predicted  value  Is  very  close  to  the  experimental  one  It  Is  still 
thought  that  this  procedure  Is  far  from  perfect,  mainly  because  the  hot 
wire  Is  soft  soldered  near  the  tip  of  the  needle  and  It  Is  this  local 
temperature  which  should  be  used.  For  each  series  of  runs  the  same 
wire  was  used,  so  that  the  error  due  to  the  support  temperature  appears 
systematically  In  all  the  data. 
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2.  Heat  Losses  and  Recovery  Temperature  as  Functions  of 

The  empirical  relation  Nu0(ReQ,M)  mentioned  In 
Is  of  the  form: 


where 


with: 


and 


with 


Nu0(Re0.M)  B  Nu0(Re0,»)$(Re0,M) 
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M,  Re  and  Kn 
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Section  3.23 


As  stated  In  Section  3.2.3,  the  above  relation  is  a  curve  fit 
for  existing  experimental  results.  This  curve  fit  gives  good  accuracy 


(±1%)  for  M  >  0.45;  below  this  value  a  more  accurate  curve  fit  Is  pro¬ 
posed  (see  (20)) 

ni  R«  -  6  9 

4>3  =  b]  +  c]  {1  -  exp[-  aj(  ReQ  -  6  )  (  Re°  -  g  )r 


where 


b1  =  0.05  +  2.9  M1 
c]  =  0.05  +  3.08  M1 


e1  =  0.6  +  1.25  M 


1.137  Jin 


e.  -  bi 

-L-l  ] 


in(l  ♦ 


el  ‘  bl 


n,  el  ‘  bi 

a]  =  0.1695  1  in[l  ♦  /  - - -  ] 


This  value  has  been  used  in  our  computation  when  M  <  0.45. 

The  relations  giving  n* ,  and  nc  are  the  following: 

1.193 


0.4930  +  Kn  '• 


nc  =  1  -  0.05  ( 


1.175  ♦  ' 


nf  -  nc  ■  0.2167  ( 


0.8521  +  M‘ 


As  stated  in  Section  3.23,  the  relations  Nu0(Re0*M)  could  be  transformed 
to  Nu  (M,T«.,p)  because  we  can  write 
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and,  substituting  for  y,R  and  using  Sutherland's  viscosity  law,  we 
get: 
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